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Abstract 

GdFes(B03)4 exhibits a structural phase transition at 156 K, antiferromagnetic order of the 
Fe 3+ moments at 36 K followed by a spin reorientation phase transition at 9 K. The reorientation 
phase transition is studied through dielectric, magnetic and heat capacity measurements under 
the application of external magnetic fields of up to 25 kOe. The dielectric constant indicates the 
existence of two distinct anomalies at T$r = 9 K that separate in temperature under external 
magnetic fields. The spin rotation phase transition is proven to be of the first order nature through 
the magnetic analogue of the Clausius-Clapeyron equation. Magneto-dielectric effect of up to 1 
percent is observed at 8 K and 7 kOe. Uniaxial magneto caloric effect along the c-axis is observed 
below the spin reorientation phase transition of 9 K. 

PACS numbers: 75.30.Kz, 75.50.Ee, 75.80.+q, 77.80.-e 
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GdFe 3 (B03)4 belongs to the trigonal system with space group R32. It is similar to huntite 
CaMg 3 (003)4, a trigonal trapezohedral structure that is one of the five trigonal types. Bo- 
rate crystals of this category are appealing because of their possible applications as single 
crystal minilasers due to their good luminescent and nonlinear optical properties. In particu- 
lar for GdFe 3 (B0 3 ) 4 single crystals, recent studies have focused on the better understanding 
of its optical properties and at the same time on its magnetic properties through phase 
matching of absorption and second harmonic generation spectra.^ The magnetic properties 
of the rare-earth iron borates are also of fundamental interest because of the existence and 
mutual interference of two magnetic subsystems (Fe and Gd). It is of fundamental impor- 
tance the understanding of the magnetic orders of the gadolinium and iron sublattices and 
the coupling between the iron spins and the gadolinium moments which contributes to the 
crystal's rich magnetic properties. 

The crystal structure of huntite has been analyzed and described elsewhere 3 -^, it consists 
of GdOe bipyramids and Fe06 octahedrons. The octahedrons form 3-fold helicoidal chains 
along the c-axis. The bipyramids are located between nearly three equal distant octahedrons. 
Rare-earth iron borates, RFe 3 (B0 3 )4 (R=Eu to Ho and Y), undergo a structural transition 
at higher temperatures and an antiferromagnetic (AFM) transition involving the Fe spins at 
TAr?a35 K.— For R=Gd a weakly first order structural phase transition at Ti=156 K changes 
the structural symmetry from R32 to P3i2!. 5 A second order phase transition at Tjy=36 
K results in the AFM ordering of the Fe 3+ magnetic moments aligned in the basal plane. 
At Ts^=9 K, another sharp phase transition occurs characterized by a spin reorientation of 
the Fe 3+ magnetic moments by 90 degrees from the basal plane to the c-axis. 

The coupling between the Gd moments and the Fe spins at low temperatures is strong in 
GdFe 3 (B0 3 )4 and it was speculated that the reorientation of the AFM iron spin system is 
triggered by the anisotropy of the Gd moments aligned with the c-axis.— This exchange in- 
teraction is indirect and involves other ions such as oxygen in the structure. Below Tjy, with 
decreasing temperature, the magnetic order grows and that causes strain via the spin-lattice 
interaction. At the spin reorientation phase transition the magnetic order experiences a sud- 
den change which is reflected in distinct anomalies of the dielectric constant, e, as recently 
discussed in several rare-earth manganites. 7 '- The dielectric constant can be measured with 
extraordinary precision and it monitors subtle changes of the magnetic system. We have 
therefore searched for dielectric anomalies at the magnetic phase transitions in GdFe 3 (B0 3 )4 
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and its correlation with magnetic and heat capacity data. Below TV, close to Tsr, we found 
two distinct anomalies of e that separate in temperature with applied magnetic fields. A 
sharp drop of e occurs at the reorientation of the Fe spins whereas a maximum of e(T) above 
Tsr indicates the onset of magnetic field induced ferroelectric order. The spin reorientation 
phase transition is proven to be of the first order nature through the magnetic Clausius- 
Clapeyron equation. Lastly, a phase diagram is proposed for magnetic fields parallel or 
perpendicular to the crystalline c-axis. 

The single crystal GdFe3(BC>3)4 was grown as described elsewhere.— The transparent crys- 
tal (shining green in color) was analyzed and oriented using a GADDS x-ray diffractometer. 
The heat capacity data was acquired via Quantum Design's Physical Property Measurement 
System (PPMS) under the application of magnetic fields of up to 6 kOe both along a- and 
c-axis (we use the hexagonal coordinate system in which ale). DC magnetization data was 
measured employing Quantum Design's Magnetic Property Measurement System (MPMS) 
in magnetic fields parallel to a and c. Smaller portions of the crystal were cut from an 
original crystal in order to align the two desired crystallographic orientations for dielectric 
investigations. Silver paint was used as electrodes for the dielectric measurements and the 
sample was mounted onto a home-made capacitance probe that was adapted to the PPMS. 
The capacitance was measured by the high precision capacitance bridge AH2500A (Andeen 
Hagerling) operating at a frequency of 1 kHz. The temperature and magnetic fields were 
controlled by the PPMS when the dielectric constant was measured. 

The three known phase transitions are reflected in different anomalies of the dielectric 
constant (Fig. 1). At Ti=156 K, the dielectric constant drops significantly since the tran- 
sition is a structural one. Ti does not depend on magnetic field. The magnetic transition 
into the AFM2 phase at Tjv is distinguished by a change in slope of e c and a minimum of 
e a at 36 K. Below T^r, e a starts increasing with decreasing T and develops the peak-like 
feature at Tsr, the transition into the AFM1 phase, as shown in Fig. 1. This is a clear 
signature of lattice softening and it has to be correlated with the changes in the magnetic 
order at and below Tjv. A similar, but positive slope change at Tjv is observed in e c . This 
behavior is analogous to the magnetic data which will be discussed later. In a system where 
there are two magnetic subsystems present with one of them being a d-metal and the other 
one an f-metal subsystem, the f-metal is magnetically polarized upon the ordering of the 
d-metal subsystem.^ In GdFe3(BOs)4 the gadolinium develops an AFM order along the c- 
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axis through inter planar exchange interactions mediated by the Gd-O-Fe bonds below Tjy. 
The dielectric anomalies observed at the magnetic phase transitions are an indication of 
strong spin-lattice interactions and are attributed to magneto-elastic effects similar to those 
observed in some rare earth manganitesii and a strong coupling of the dielectric response 
with the rare earth magnetic moment was also reported in these compounds. 7 ' 8 The softness 
of the lattice (as expressed by the sharp increase of e a (T)) is therefore an indication of the 
softening of the magnetic order as T approaches Tsr resulting in the Fe-spin reorientation 
triggered by the Gd moments and the abrupt decrease of e a . Thermal hysteresis was ob- 
served in e(T) above Tsr at zero magnetic field. The hysteretic behavior of e (T) extends 
over several degrees just above Tsr- Upon cooling in zero field the dielectric constant below 
T N increases continuously for T — ► T S r and drops rapidly right below Tsr- However, upon 
heating, after the sharp increase at Tsr there appears a distinct maximum of e a at Tm about 
0.7 K higher in temperature (Fig. 2). This second anomaly at T^ needs to be explored in 
more detail. 

The external magnetic field, H, shows an interesting effect on both dielectric anomalies 
(Tsr and Tm)- If H is aligned with the c-axis the spin reorientation transition is suppressed 
and Tsr quickly decreases with H c in accordance with recent magnetic data.— However, T^ 
is barely affected by H c and the separation of both anomalies increases in magnetic fields 
as shown in Fig. 2. In contrast, for magnetic fields parallel to a, Tsr remains constant but 
Tm shifts to higher temperature (Fig. 3). At the same time the enhancement of e a (T) is 
largely reduced and disappears for fields exceeding 7 kOe. This qualitative behavior of the 
dielectric constant in magnetic fields is closely correlated with the magnetic properties of 
the coupled system of the Fe spins and the Gd moments. It is interesting to note that recent 
measurements of the magnetic field-induced electric polarization, P, of GdFes(B03)4 have 
shown a sudden increase of P in crossing T^ with increasing field.— This may indicate the 
existence of a field-induced ferroelectric phase (FIP) and it was proposed that higher order 
magneto-electric couplings are responsible for the observed phenomena. From our dielectric 
measurements we can uniquely identify the boundaries of this new phase by the sudden 
drop of e a (T) at T S r and the distinct maximum of e tt (T) at the high temperature end (T M )- 
Therefore, the phase diagram of GdFe3(B03)4 includes three distinct phases (AFM2, AFM1, 
and FIP) below the AFM ordering temperature, Tjv, of the iron system. 

The results of magnetization measurements with fields oriented parallel to a and c are 
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summarized in Fig. 4. Surprisingly, the a- and c-axis DC susceptibilities are equal above 
TV within the experimental resolution, i.e. GdFe3 (1303)4 is magnetically isotropic at high 
temperatures although there is evidence for easy plane anisotropy of the Fe spins and easy 
axis anisotropy of the Gd spins.— The anisotropies of both magnetic ions are obviously 
correlated resulting in a complete isotropic balancing just like some rare-earth metal ferrite- 
garnets and orthoferrites.— The high temperature Curie-Weiss extrapolation yields a Curie- 
Weiss temperature of Qcw = — 34 K and an effective magnetic moment fi e ff = 12.2 fis in 
good agreement with other reports^ and the theoretically expected value of (JL e ff = 12.96 
(Gd 3+ , S=7/2 and Fe 3+ , S=5/2). Deviations from the Curie- Weiss behavior due to short 
range AFM correlations are detected below 100 K. 

Below Ttv the magnetic response becomes anisotropic and depends on the direction of 
the probing field. For H || c the DC susceptibility, x c (T), does not exhibit any anomaly at 
Tjv and it is independent of the magnetic field (up to 10 kOe) between T N and T S r- Xc 
decreases suddenly at Tsr as a result of the spin reorientation aligning the Fe-spins with 
the c-axis (Fig. 4). The major effect of H c is the shift of Tsr to lower T and eventually the 
suppression of the spin reorientation for fields above 8 kOe. It is obvious that the magnetic 
response strongly depends on the angle between the external field and the Fe spins. This is 
justifiable since it is energetically unfavorable to keep the Fe spins aligned with the external 
magnetic field. This qualitative behavior is reflected in the field dependence of the dielectric 
constant (Fig. 2). e (T) increases below T/v and exhibits a distinct maximum at Tm~10 
K. The sharp drop of e a (T) at lower temperature coincides with Tstj(H) as determined 
magnetically. e a (T) is only weakly dependent on the field H c between T^ and Tm as is the 
magnetic susceptibility, Xc- Whereas the value of T M is not affected by the field oriented 
along the c-axis Tsr(H) decreases rapidly reaching zero at about 8 kOe and separating 
the two dielectric anomalies in the magnetic field. The peak of e a (T) indicating the phase 
boundary between the AFM2 phase and the FIP phase exists well above 8 kOe, the critical 
field that suppresses the AFM1 phase. Details of the phase diagram will be discussed later. 

The magnetic response to an in-plane magnetic field, H a , is completely different. At 
low magnetic fields the susceptibility Xa levels off below T^v and is clearly lower than Xc 
(Fig. 4). This can be attributed to the AFM order of the Fe spins below T^ reducing 
the susceptibility for H || a. Close to T S r, however, \a increases again and its temperature 
dependence below Tsr. appears to be a continuation of the c-axis susceptibility, Xc (discussed 
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above). This behavior agrees well with the observation that the magnetic susceptibility is 
little affected by the AFM order if the field is perpendicular to the Fe spins. At T$r the Fe 
spins rotate towards the c-axis and align themselves with the magnetic field. Xa(H) rapidly 
increases with H in the temperature range Ts^<T<IV and it approaches the values of Xc 
above 7 kOe. Tsr is independent of H a , as is also reflected in the dielectric data of Fig. 3 
(the sharp drop of e a (T) marks the spin reorientation transition). e a (T) between Tsr and 
T^, however, is strongly affected by H a and its enhancement below is reduced by the 
magnetic field. Thereby the maximum of e a (T) is shifted to higher temperatures (Fig. 3) 
reaching the Neel temperature at a relatively low field of 4 to 5 kOe. This limits the AFM2 
phase and extends the temperature range of the FIP phase to Ttv=36 K. 

The dielectric anomalies observed in GdFe3(BOs)4 at the magnetic transitions reveal an 
interesting correlation between the different magnetic subsystems (Fe and Gd) and the spin- 
lattice coupling leading to the modifications of e(T). To arrive at a deeper understanding 
the magnetic correlations have to be considered in more detail. The phase diagram and the 
magnetic order of GdFe3(BOs)4 was recently investigated by AFM resonance experiments 
and it was suggested that the Fe spins at high temperatures are aligned in the a-b plane 
(easy plane anisotropy) whereas the Gd moment experience a strong easy axis anisotropy 
orienting them along the c-axis.— The AFM order of the iron moments below T^v couples to 
the gadolinium moments resulting in an AFM alignment of the Gd spins ("magnetic polar- 
ization"). With decreasing temperature the Fe-Gd coupling grows stronger and eventually 
triggers the reorientation of the Fe spins resulting in a transition from the non-collinear 
(Ts#<T<Tyv) to a collinear (T<Tst,;) magnetic order of both subsystems along the c-axis. 
This phenomenon, related to the strong single ion anisotropy of rare earth ions, is not un- 
common and it was observed in a number of different compounds involving rare earth ions 
and d-elements . 15 i 16 i 17 Depending on the strength of the anisotropy factors the coupling be- 
tween the f-moments and d-spins results in the rotation of the d-spins towards the f-moments 
(such as in R2CUO4, R=Pr, Nd, Sm, Eu) 1 ^ 1 ^ or vice versa (e.g. in orthorhombic HoMnOs).^ 
In rare cases the f-d exchange can result in an in-plane rotation of the d-spins as observed 
in hexagonal HoMnOs.— Dielectric anomalies associated with spin reorientation transitions 
have been observed in several manganites^S and the significance of the spin-lattice coupling 
and magnetoelastic effects have been shown.— The coupling between the f- and d-spins 
forming a 90 degree angle can be mediated by the antisymmetric Dzyaloshinskii-Moriya 



6 



interaction (that is proportional to the vector product of f- and d-moments and gives rise 
to weak ferromagnetism) or by a pseudodipolar interaction arising from the anisotropy of 
the f-d exchange.— In any case, the exchange interaction maximizes if the two moments are 
perpendicular to each other. 

The Fe spin reorientation observed in GdFe 3 (B0 3 )4 at Tsr is triggered by this f-d ex- 
change mechanism. We propose that the observed enhancement of the dielectric constant 
below TV (Figs. 1 to 3) is closely related to the indirect f-d exchange coupling and the 
magnetic polarization of the Gd moments. The increase of e a (T) as T approaches Tsr at 
H=0 reflects the increased coupling strength and AFM Gd sublattice polarization. The 
spin-phonon coupling causes the lattice to soften as the magnetic order becomes unstable 
for T — ► T S r- The combined effects on e a saturate at T 5R (H=0) and it disappears abruptly 
with the spin reorientation of the iron at lower T resulting in the sudden decrease of the 
dielectric constant. With the magnetic field applied along the c-axis Tsr is lowered but the 
polarizing effect on the Gd moments and the enhancement of e a (T) still reaches its max- 
imum at 9 K explaining the peak of e a remaining H independent at Tm~9 K. With the 
magnetic field applied along the a-axis the Gd moments are tilted away from the c-axis and 
the in-plane AFM order of the Fe-spins is gradually reduced, as can be concluded from the 
field effect on the susceptibility, \a (Fig. 4). Both effects reduce the f-d exchange coupling 
and the magnetic polarization of the Gd-spins and its effect on the lattice is reduced. The 
saturation of the lattice softness happens at higher temperature in the in-plane magnetic 
field, the maximum of e (T) is suppressed, and Tm increases with H (Fig. 3). Although 
this discussion is qualitative it does explain the two anomalies observed in the temperature 
dependence of the dielectric constant and their opposite dependence on the magnetic fields 
parallel or perpendicular to the c-axis. The softening of the lattice at Tm eventually re- 
sults in a ferroelectric displacement as observed in recent polarization measurements^ that 
indicate the existence of a sizable field-induced polarization between Tsr and Tm for both 
orientations (c-axis and in-plane) of the magnetic field. The corresponding magnetostriction 
data provide the experimental prove for the lattice anomalies at the phase boundaries of the 
FIP phase, T S r and Tm- 

The thermodynamic signature of the magnetic phase transitions at and Tsr is given 
by sharp peaks of the zero-field heat capacity^ at the phase transitions. In external magnetic 
fields, H || c, the heat capacity peak at Tsr shifts to lower T , in accordance with the magnetic 
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and dielectric data discussed above. Fig. 5 shows the heat capacity close to Tsr at different 
values of H c . All peaks are relatively sharp and symmetrical indicative of a first order phase 
transition. The heat capacity peak does not shift with H || a which is consistent with our 
dielectric and magnetic data. From the structural point of view the exchange interaction 
between gadolinium ions should be weak since they are interconnected by BO3 triangles that 
are 2.4 A apart. Thus Tm is not observed in the heat capacity or magnetization data but 
only in dielectric data. 

The magnetic phase diagram of GdFe3(BC>3)4 at lower temperatures derived from the heat 
capacity, magnetization and dielectric measurements is shown in Fig. 6 for both orientations 
of the magnetic field. Tsr is not affected by H a , however, it decreases quadratically with H c , 
in agreement with the results of recent resonance^ and magnetoelectric experiments. 12 The 
field dependence of Tm is also shown in Fig. 6 by dashed lines. Tm merges with T^ for H a « 
5 kOe in good agreement with the anomalies detected in polarization and magnetostriction 
measurements. However, for H || c, Tm(H) is almost constant and extends to far higher 
fields than Tsr- Tsr and Tm form the phase boundaries of the FIP phase where magnetic 
field-induced polarization was observed.— It is interesting to note that the heat capacity 
does not show any anomaly at Tm as it does at Tsr- The transition across Tm is obviously 
not accompanied by anomalies such as a sizable change of volume etc. This indicates a more 
subtle change in the correlated magnetic systems of Gd and Fe spins coupled to the lattice 
that results in the distinct peak of the dielectric constant. 

Given the sharp nature of the heat capacity peak anomaly at Tsr as well as its sudden 
drop in magnetization and dielectric constant, a first order phase transition is presumed to 
occur. To further verify the validity of the first order nature at Tsr we make use of the 
magnetic analog of the Clausius-Clapeyron equation, 

-AS(B)/AM(B) = dB/dT 5i? (B) 

where AS is the change in entropy, AM is the change in magnetization and dB/dTsR 
is the inverse slope of Ts#(B) with respect to the internal magnetic field B. The set of 
values for AS are obtained from integrating the area underneath the peak anomaly of Cp/T 
at different fields (Fig. 5). The set of values for AM is obtained from the change of the 
magnetization at Tsr (Fig. 4). dB/dT sr is extracted from the phase diagram of Fig. 6 and 
B=/i (H+M). For the case of H || c, equation (1) is well fulfilled as shown in Fig. 7. This 
proves the first order nature of the spin reorientation transition. 
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In conclusion, we have demonstrated a strong correlation between the magnetic order 
and the dielectric properties of GdFe 3 (B0 3 )4. Below the AFM ordering temperature of the 
Fe-spins, T^, e a (T) increases and reaches a maximum at the spin reorientation transition 
temperature, Tsr, in zero magnetic field. We interpret the apparent lattice softness as 
a consequence of the indirect exchange coupling between the Gd- and Fe-spins. For non- 
zero magnetic fields the e a (T)-maximum (at Tm) and the spin reorientation transition are 
separated and appear at different temperatures. The spin reorientation transition is proven 
to be a first order phase transition through the magnetic analogue of the Clausius-Clapeyron 
equation. 
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FIG. 1: (Color Online) Dielectric data showing all three major anomalies. Ti=156 K, TV=36 K 
and Tsr=9 K. Inset shows the behavior of e(T) around Tat when e is measured along the hexagonal 
a- and c-axis for H=0. 

FIG. 2: (Color Online) Splitting of the spin reorientation phase transition under fields parallel to 
the c-axis. Tsr shifts lower while Tm stays the same. 

FIG. 3: (Color Online) Splitting of the spin reorientation phase transition under fields parallel to 
the a-axis. Tsr remains constant while Tm shifts higher. 

FIG. 4: (Color Online) Mass susceptibility data at different magnetic fields parallel and perpen- 
dicular to the c-axis. Tsr shifts lower and AM increases under magnetic fields parallel to the 
c-axis (open). Tsr stays the same if magnetic fields are aligned with the a-axis (solid). Xa is only 
affected between Tsr<T<Tn. 

FIG. 5: (Color Online) Heat capacity data at different magnetic fields parallel to c-axis. Tsr shifts 
lower and AS decreases with H c . 

FIG. 6: (Color Online) Phase diagram of GdFe3 (603)4 at low temperatures including Tsr, Tm and 
Ttv derived from dielectric data. AFM1 is the state where the Fe spins are ordered antiferromagnet- 
ically along the basal plane, AFM2 is the state where the Fe spins are ordered antiferromagnetically 
along the c direction. FIP is the state where there is field induced polarization. 

FIG. 7: (Color Online) The Clausius-Clapeyron equation satisfied at different magnetic fields H c . 



11 




H — > — i — ' — i — ' — i — > — i — ' — i — ' — |— 
5 10 15 20 25 30 

T(K) 



30- 



20- 



AFM2 



▲ 


T 




N 




T 




M 


— ■ 


T 




SR 




AFM1 



"T 

8 







2 4 6 

H (kOe) 



10 





— 1 1 1 1 1 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 

B (Tesla) 



